Excess dissipation in a single-electron box: The Sisyphus resistance 
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We present measurements of the ac response of a single-electron box (SEB). We apply an rf signal 
with a frequency larger than the tunneling rate and drive the system out of equilibrium. We observe 
much more dissipation in the SEB then one would expect from a simple circuit model. We can 
explain this in terms of a mechanism that we call the Sisyphus resistance. The Sisyphus resistance 
has a strong gate dependence which can be used for electrometery applications. 



Dissipation in quantum systems has been an impor- 
tant area of research for many years. It has lately gained 
renewed attention as quantum systems have increasingly 
come to the forefront of technology, including research 
in nanotechnology and quantum information. In many 
of these contexts, dissipation can be modeled using a 
simple two-level system (TLS) or an ensemble of them. 
It is therefore generally useful to study dissipation in 
TLSs. For instance, in the context of quantum infor- 
mation where any dissipation causes unwanted de- 
coherence, TLSs are important in several different ways. 
First, the basic building block of all quantum bits is an 
effective TLS. In addition, the presence of parasitic TLSs 
in dielectrics has been shown to result in losses which de- 
grade the quality factor of electrical resonators and limit 
the coherence of superconducting qubits [2] . Fast driving 
of a TLS through an avoided level crossings (ALC) re- 
cently received considerable attention, and has been an- 
alyzed in terms of Landau-Zener (LZ) transitions 
and dressed states 0. Recently, in a setup related to the 
one studied here, Sisyphus cooling [7] and amplification 
was observed in a superconducting circuit [8] . 

In this Letter, we study a mesoscopic circuit consisting 
of a small metallic island connected by a tunnel junction 
to a much larger reservoir. Charging effects (Coulomb 
blockade) result in well defined energy levels which be- 
come degenerate at a specific bias point. An ac drive 
is used to cyclically drive the system through this level 
crossing. Due to the low transparency of the tunnel bar- 
rier, the coupling between the levels is negligible and the 
probability of a LZ transition when crossing the degen- 
eracy point is very close to unity. However, due to the 
large degeneracy of the electronic states on the island, the 
total system can have a significant relaxation (or tunnel- 
ing) rate. If the frequency of the drive is comparable to 
the relaxation rate of the system, alternating excitation 
and relaxation of the system lead to excess dissipation 
which can be directly measured. We call this process the 
Sisyphus resistance. We develop a quantitative model of 
the behavior that shows very good agreement with the 
measured response. 
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The charging effects which give rise to the quantization 
of the energy levels have a very peculiar effect on dissipa- 
tion in the circuit. Far away from the degeneracy point, 
the Coulomb blockade prevents tunneling, and results in 
a high Sisyphus resistance and low dissipation. How- 
ever, at the degeneracy point, the charging effects result 
in a low Sisyphus resistance with dissipation much larger 
than that expected from a gate capacitance in series with 
an Ohmic resistor with the value of the tunneling resis- 
tance. Fig. \Vp. We note that the Sisyphus mechanism 
described here is not limited to Coulomb blockade de- 
vices, but should be a generic property of any system 
with an energy level crossing. In an analogous way to 
the Sisyphus mechanism discussed here, it was suggested 
that phonon pumping of TLSs can dominate losses in 
micro-mechanical resonators at low temperatures Q and 
one could imagine having the same effect in electrical res- 
onators. It has been suggested that TLSs may dominate 
the loss in electrical resonators under conditions impor- 
tant for quantum information applications [2|. 

The circuit of interest in this Letter is very similar to 
the radio-frequency (rf) SET ^ except that the SET 
is replaced by a single- electron box (SEB) (see Fig. [1^). 
The SEB consists of an aluminum island connected by a 
tunnel junction to a charge reservoir and three different 
capacitive gates. There is an rf gate which couples the 
SEB to an on-chip lumped-element resonator. The SEB 
also has a dc gate to adjust the static potential of the 
island and a microwave (mw) gate used for spectroscopic 
characterization of the SEB. We perform rf reflectometry 
measurements by sending a continuous rf signal to the de- 
vice. On resonance, the rf signal will excite the resonator, 
which in turn will drive the potential of the SEB island 
around the working point set by the dc gate. Depending 
on the dissipation in the resonator caused by the SEB, 
the magnitude of the reflected signal will change. We 
measure both the magnitude and phase of the reflected 
signal to characterize this dissipation. 

To explain the observed dissipation, we can consider 
the two energy levels, Eq and Ei, in Fig. [TJl. The energy 
of each level can be controlled by the parameter Ug, the 
normalized gate voltage in units of induced electrons. At 
the degeneracy point, Ug = 0.5, the two levels cross. Let 
us now assume that the system is biased at a point nPg and 
driven around this point by a fast rf drive of amplitude 
Sug. We imagine following the system through one cycle 
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FIG. 1: (a) A scanning-electron micrograph of the SEB. The 
island of the SEB is 8 fim long and 100 nm wide. The SEB 
has three different gates: one rf gate connecting the SEB to 
a parallel resonator used for readout, one dc gate to set the 
working point and one microwave gate used in the character- 
ization. To the right of the micrograph are different circuit 
models used in the discussion of the measured response, (b) 
The circuit model used in analyzing the response. The SEB 
is modeled by the voltage dependent resistor Rsis- The in- 
ternal losses in the resonator, mostly due to the normal top 
layer, are modeled as a parallel resistor Rioss- The coupling 
capacitor Cc is used to decouple the resonator from the 50 Q. 
impedance of the environment and thereby increase the qual- 
ity factor of the resonator, (c) A naive model of the SEB, 
which is compared with the Sisyphus resistance in the text. 
Cj is the geometric capacitance of the tunnel junction, (d) 
The process behind the Sisyphus resistance. -Bo and Ei are 
the electrostatic energies of the two charge states. Starting 
from a charge state of lower energy, the rf signal Sug drives the 
system past the degeneracy point and into an excited state, at 
the same time charging the system. When the system eventu- 
ally relaxes by tunneling, the charging energy is lost, (e) The 
tunnel rates between the two charge states as a function of the 
normalized gate voltage, rig, for three different temperatures. 



of the rf drive (see Fig. [TJl) . In the first half cycle, we 
start on the left (right) side of the degeneracy point in 
the ground state. Moving across the degeneracy point 
at a high rate the system stays in the same state, which 
becomes the excited state. As we move away from the 
degeneracy point, the tunneling rate, r+^.-j, for electrons 
on to (off from) the island, increases (see Fig. [1^) until 
a tunneling event occurs. The energy that had been put 
into the system is now dissipated and we have to start 
over again in the second half cycle. In contrast, when we 
are far from degeneracy the energy put into the system 
during the first half cycle is given back during the second. 

The electrostatic energy of the SEB is £^ = Ec{n—ng)^ 
where n is the number of electrons added to the is- 



land from its neutral state and Ug is the charge induced 
on the gate capacitance by the external voltage. Here 
Ec = e^/SCs is the charging energy and Cs is the total 
capacitance of the SEB island. For low enough temper- 
ature, ksT <^ Ec, and a gate charge ranging between 
and 1, only two charge states {Eq and Ei) need to 
be considered. The dynamics of the driven SEB is then 
governed by the following master equation (ME) 

Po = r_Pi - r+Po (la) 
Pi - r+Po - r_Pi (lb) 

where Pq and Pi are the probabilities of being in the two 
charge states and the tunneling rates are given by the 
orthodox theory for single-electron tunneling [ll| 

TAE/h Rk 

and where AE =^ Ei - Eq ^ Ec (1 - 2ng{t)), ng{t) = 
n^g -f Sug sm{uj{)t) and 27r/o is the angular frequency 
of the rf drive. Here ?i° — Cg'^Vg/e is the dc gate charge, 
Px = h/e^ is the resistance quantum and Rt the tunnel- 
ing resistance of the junction in the SEB. The normalized 
rf amplitude is Sug = Cg^V^^ /e, where Vg^ is the volt- 
age amplitude inside the resonator. In order to solve the 
ME, we expand the tunneling rates ([2]) around the work- 
ing point, n^g, to first order in 5ng and insert them into 
the ME. We then linearize the equation to get a first or- 
der linear differential equation for Pi (or Pq) that can be 
solved analytically to get Po(i) and Pi{t). The average 
power dissipation, i.e. the energy transfered to the bath 
per cycle of the rf drive, T = l//o, can then be calculated 
as 

Psrs f dt [P^T^AE - PoF+AP] (3) 
^ Jo 

For a resistor, R, driven by an AC-voltage, SV, the av- 
erage power dissipation is P = SV^ /2R. By comparing 
these two values we can define the Sisyphus resistance of 
the SEB, 

(4) 

where AE^ — Ec (l — 2ng) is the energy splitting at the 
working point and 7 = r+ -I- r_ = coth (^f^) 

is the equilibration rate of the system. This expression 
is only valid for Sug < kBT/Ec where the linear expan- 
sion of the tunneling rates is a good approximation. In 
this limit of small amplitudes, the Sisyphus resistance 
is independent of the actual amplitude. In addition to 
the Sisyphus resistance, the impedance of the SEB can 
also have a reactive component. However, our theoret- 
ical estimates predict that the gate dependence of this 
reactance is negligible in our experiment. 
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The device was fabricated in a multilayer process. 
Starting from a high- resistivity silicon wafer with a native 
oxide, the wafer was first cleaned using rf back sputter- 
ing directly after which a 60 nm thick layer of niobium 
was sputtered. To pattern the niobium, we used an Al 
mask made by e-beam lithography and e-beam evapo- 
ration. The niobium was then etched in a CF4 plasma 
to form the inductor and bottom plates of the capaci- 
tors (see Fig. [1^). The Al mask was removed with a 
wet-etch based on phosphoric acid. Using PE-CVD, we 
then deposited an insulating layer of 200 nm of silicon 
nitride which covers the whole wafer; connections to the 
niobium layer were only made through capacitors. After 
using a combination of DUV photolithography to define 
bonding pads along with e-beam lithography to define 
the top layer of the capacitors, a 3/80/10 nm thick layer 
of Ti/Au/Pd was deposited by e-beam evaporation. Fi- 
nally, the layer containing the SEB was made by e-beam 
lithography and two-angic shadow evaporation of 30-1-50 
nm of aluminum, with 6 min of oxidation at 4 mbar. 

The device was cooled in a dilution refrigerator with 
a base temperature of about 20 mK. For the readout, 
we used an Aeroflex 3020 signal generator to produce 
the rf signal. The signal was heavily attenuated and fil- 
tered and was fed to the tank circuit via two Pamtech 
circulators positioned at the mixing chamber. The re- 
flected signal was amplified by a Quinstar amplifier at 4K 
with a nominal noise temperature of 1 K. The in-phase 
and quadrature component of the signal were finally mea- 
sured using an Aeroflex 3030 vector digitizer. 

In order to determine various device parameters, we 
performed dressed-state spectroscopy of the box in the 
superconducting state T he q uantum capacitance of 
the box was used for readout In order to deter- 

mine the charging energy, we applied microwave frequen- 
cies of =10-11 GHz to the mw gate while slowly sweep- 
ing the dc-gate. Multiphoton resonances occur when 
"nihff^ 4i?c'(l — Ug). From the positions of the multi- 
hoton resonances, we can determine the charging energy 
For this device, Ec/h — 15.0 ± 0.1 GHz which cor- 
responds to a total capacitance of Cs — 1.29 fF for the 
SEB. Using this value of Ec, we then extracted a value 
for the Josephson coupling energy of Ej/h = 3.6 ± 0.2 
GHz from conventional spectroscopy. 

For the normal state measurements, we applied a par- 
allel magnetic field of 600 mT in order to suppress su- 
perconductivity in the Al but not in the Nb. The res- 
onator circuit was first characterized by measuring the 
reflection coefficient, 5*11, as a function of frequency. Fit- 
ting the real and imaginary part of Sn simultaneously, 
the parameters of the resonator could be determined. 
From the fit we extracted L = 322 nH, Cc= 92.1 fF, 
C=109.5fF and i?;oss=608kri. This corresponds to a 
resonant frequency of /o=624.7MHz and a total Q- value 
of 97. The total Q-value has contributions both from 
the coupling to the external impedance of i?o = 50 ft, 
Qq = {C + Cc)/C^(jJoRo = 121, and from internal losses, 
Qr = ujoR{C+Cc), where R = RiossRsis/{Rioss + Rsis)- 
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FIG. 2: The reflected amplitude, for low rf amplitudes (Sug « 
0.04), as a function of the the bias point, n^, for three differ- 
ent temperatures. The theory (solid lines) are fit to the mea- 
sured data (dots) using the analytic formula for the Sisyphus 
resistance, Q, inserted into the expression for the reflection 
coefficient of the combined system. 



Away from the degeneracy point, where R — Riossi we 
find Qn — 481. The main contribution to the internal 
loss, which we represent by Rioss, is probably due to the 
top layer of the capacitors which was made out of gold. 
For similar circuits with superconducting Al as the top 
layer, we see much less internal loss. 

Measurements of the ac-response were performed by 
slowly sweeping the dc-gate and continuously monitor- 
ing 5*1 1 at the fixed frequency /q. There was little or no 
change in the phase of Sn but there was a substantial 
reduction in its magnitude when the SEB was close to 
its degeneracy points at Ug = ±0.5. The measurements 
were done for a range of temperatures and rf amplitudes. 
In Fig. [21 the measured response is shown for three dif- 
ferent temperatures as a function of the dc-gate charge. 
Here we used a low enough amplitude {Sug w 0.04) that 
the observed response was amplitude independent. The 
measured data were then fit to theory by replacing the 
SEB by its Sisyphus resistance using ^ and calculat- 
ing the refiection coefficient for the combined system at 
resonance. 



Sniujo) 



^l(wo) - Rq _ k {Qo - Qr) - i 



2'l(i^o) 



k{Qo + QR) + V 
1 



(5) 



where ^^(wo) is the impedance on resonance of the total 
load, as pictured in Fig. [Hd, and k = Cc/{C Cc)- The 
two fitting parameters were Rt and C'^-^, which we were 
not able to determine by independent measurements. 
The two higher temperature curves were fit simultane- 
ously, giving i?T=49kf7 and C^-^=0.31fF. The extracted 
value of Rt agrees with the measured value of Ej within 
the experimental error. The last curve was then fit by 
adjusting only the temperature, yielding an electron base 
temperature of T = 41 mK. 

We see that at the degeneracy point we get roughly 
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FIG. 3: The normalized reflected ampfltude, |5ii|, measured at the base temperature, as a function of the bias point, rig, and 
the rf amplitude, Sug. In (a) are the measured data and in (b) is the calculated response, obtained by numerically solving the 
full ME. The only adjustable parameter is the total rf coupling, which sets the y-scale. In the inset, we show three line cuts 
comparing the measured (dots) and calculated response (solid lines). We see a very good agreement between theory and data. 



17% absorption. This is a factor of eight more 
than the 2% that is expected from the naive circuit 
model shown in Fig. [T};. If we compare the dis- 
sipated power due to the Sisyphus resistance, Psis 
(Eq. ([3])), to the dissipation, expected from this 
naive circuit model, Pdr, we get that Psis/Pdr — 
\{l + {Lo^RTC^f) / {{kBT/Ec? + {uJoRTC^f) ^ 8. 
A simple way to understand the large discrepancy is to 
consider the voltage that develops over the tunnel junc- 
tion in the two cases. In between tunneling events, the 
effective resistance of the blockaded junction is much 
higher than RT. Thus, the typical voltage developed 
across the blockaded junction is much larger than would 
develop across a junction that is not blockaded. The 
typical energy scale of the dissipative tunneling events is 
therefore larger. This leads to the curious result that, 
due to the Sisyphus effect, the dissipation in the SEB 
is actually significantly increased by the Coulomb block- 
ade, which is the opposite of what is typically found in 
single-electron devices. At the degeneracy point, there is 
a small discrepancy between the data and theory at the 
lowest temperatures. This could possibly be explained 
by temperature dependent renormalization effects 
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Using the parameters extracted from the measure- 
ments in Fig. [21 we can then calculate the response for 
arbitrary 5ng and T. In Fig. |31 we show the measured 
response at the base temperature as a function of 5ng 
together with the calculated response, obtained by nu- 
merically solving the full ME. The only additional fitting 
parameter is the total rf coupling, which sets the y-scale 
of the data. We see a very good agreement between ex- 
periment and theory. 

Since the Sisyphus resistance has such a strong gate 
dependence, it is possible to use it as the basis for a 



very sensitive electrometer. If the SEB is biased on the 
side of the degeneracy point where the response has the 
maximum slope, a small variation in the gate charge 
will give a large change in the magnitude of the re- 
flection coefficient. We have estimated the sensitivity, 
6Q = e{5VlV) \dSiildng\~' where Fand 5V are the ap- 
plied signal and noise voltage at the input of the ampli- 
fier, respectively. We calculate 5V by assuming a system 
noise temperature of 1.5 K due to noise added by the 
amplifiers and the insertion loss between the sample and 
the first amplifier. We used the calculated response in 
Fig.[3]to calculate the transfer coefficient, 5\Sii\/ Sug, as 
a function of rf amplitude. We get a best sensitivity of 74 
/j,e/\/Hz. This is not to far from the best measured sen- 
sitivity which we obtained: 86 /ie/-\/Hz. In order to op- 
timize the sensitivity, as a first approximation, we maxi- 
mize the transfer coefficient \dSii/ dng\ using Eq. ([5]) and 
the analytical expression (|4|). If we use an applied voltage 
corresponding to 5ng ~ ksT / Ec we get that 



V2 



V dx 



(6) 



In- 



where i?(n°) = Rioss\\Rs%s{n^g) and a;(n^; - ^.^j, 
serting the parameters for this sample results in a best 
sensitivity of 117/xe/\/Hz. Although there is a deviation 
from the numerical calculation, which can be attributed 
to a too low estimate of the optimal probe amplitude, 
the formula still give us valuable information on how to 
improve the sensitivity. By reducing the losses in the res- 
onator to reach an internal Q- value of 10* (which we have 
obtained in similar resonator with a superconducting top 
layer) and by decreasing the coupling capacitance, Cc, 
and the junction capacitance, Cj, by a factor three each 
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we estimate that the sensitivity can be improved by an 
order of magnitude. 

The use of an rf-SEB instead of the more common rf- 
SET has some advantages. First, it is a simpler circuit, 
requiring just one tunnel junction, which is beneficial in, 
for example, molecular electronics where junctions are 
very difficult to fabricate. In fact, in Ref. [l^ the same 
mechanism of dissipation that has been discussed here 
was used to probe the potential inside nanostructures. 
The model for the response developed in this Letter could 
be used to understand and improve the response in such 
a system. The single tunnel junction also makes the rf- 
SEB insensitive to electrostatic discharges, since there 
is no dc path through the device. This makes the rf- 
SEB more suitable for applications in demanding envi- 



ronments, such as for scanning probes fl6|. 

In conclusion, we have measured the dissipation in a 
single-electron box driven by an rf signal. The observed 
dissipation is surprisingly large and cannot be explained 
by a simple circuit model. We explain this result using 
a master equation description of the Sisyphus resistance. 
We also demonstrate that this phenomenon can be used 
for electrometry. 
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